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Abstract 

The Hohenberg-Kohn (HK) theorems of bijectivity between the external scalar potential and 
the gauge invariant nondegenerate ground state density, and the consequent Euler variational 
principle for the density, are proved for arbitrary electrostatic field and the constraint of fixed 
electron number. The HK theorems are generalized for spinless electrons to the added presence 
of an external uniform magnetostatic held by introducing the new constraint of hxed canonical 
orbital angular momentum. Thereby a bijective relationship between the external scalar and 
vector potentials, and the gauge invariant nondegenerate ground state density and physical current 
density, is proved. A corresponding Euler variational principle in terms of these densities is also 
developed. These theorems are further generalized to electrons with spin by imposing the added 
constraint of hxed canonical orbital and spin angular momentum. The proofs differ from the 
original HK proof, and explicitly account for the many-to-one relationship between the potentials 
and the nondegenerate ground state wave function. 
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I. INTRODUCTION 


Q 


The Hohenberg-Kohn (HK) theorems [jj constitute a fundamental advance in quantum 
mechanics. As a consequence they have furthered our understanding of the electronic struc¬ 
ture of matter: atoms, molecules, solids, clusters, surfaces, lower dimensional electronic 
systems such as heterostructures, quantum dots, graphene, etc. Matter, according to HK, is 
described as a system of N electrons in an external electrostatic held ^^(r) = —Vn(r). The 
hrst HK theorem dehnes the concept of a basic variable of quantum mechanics. Knowledge 
of this gauge invariant property - the nondegenerate ground state density p(r) ~ is of two-fold 
signihcance: (a)\i determines the Schrodinger theory wave functions T of the system, both 
ground and excited state; (b) As the wave function T is now proved to be a functional of the 
basic variable, it constitutes together with the second HK theorem - the energy variational 
principle for arbitrary variations of the density - the basis of theories of electronic structure 


such as of Hohenberg-Kohn |l|, Kohn-Sham 


2j, and quantal density functional theory 


M- 


The theorems are valid for arbitrary confining potential n(r) and electron number N, but 


are derived for the constraint of fixed N. In this paper we generalize the HK theorems 
for spinless electrons to the added presence of an external uniform magnetostatic field. As 
the presence of the magnetic field constitutes a new degree of freedom, we introduce the 
further natural constraint of fixed canonical orbital angular momentum. Thereby we prove 
that the basic variables in quantum mechanics in a uniform magnetic field are the gauge 
invariant nondegenerate ground state density p(r) and physical current density j(r). These 
theorems are then further generalized to electrons with spin by imposing the constraints of 
fixed canonical orbital and spin angular momentum. 

The generalization is motivated by the considerable recent interest in yrast states which 
are states of lowest energy for fixed angular momentum. These states have been studied 
experimentally and theoretically for both bosons and fermions, e.g. rotating trapped Bose- 
Einstein condensates jb], and harmonically trapped electrons in the presence of a uniform 
perpendicular magnetic field 7j . The theorems derived are applicable to all experimentation 
with a uniform magnetic field such as the magneto-caloric effect j^, the Zeeman effect, 
cyclotron resonance, magnetoresistance, the de-Haas-van Alphen effect, the Hall effect, the 
quantum Hall effect, the Meissner effect, nuclear magnetic resonance, etc. 

The manner by which a basic variable is so defined is via the proof of the first HK 
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theorem for v-representable densities. To explain this, and to contrast the present proofs 
with the HK proof, we hrst briefly describe the HK argnments. The HK theorems are 
proved for a nondegenerate gronnd state. Particnlarizing to electrons withont any loss of 
generality, the Hamiltonian H in atomic nnits (charge of electron —e; \e\ = h = m = 1) 
is iP = I J2kPk + ^ Si e 1/kfc — ^i\ + Sfc where the terms correspond to the kinetic 

T (with momentnm p^. = —iVrj,), the electron-interaction IT, and external potential V 
operators, respectively. The Schrodinger equation is hf(R)\h(X) = ii^T(X), where 
are the eigenfunctions and eigenenergies, with R = ri,...,r 7 v;X = Xi,...,X 7 v;x = rcr being 
the spatial and spin coordinates of the electron. The energy E is the expectation E =< 
>I>(X)|//(R)|>I>(X) >. In the first HK theorem it is initially proved (Map C) that there is a 
one-to-one relationship between the external potential n(r) and the nondegenerate ground- 
state wave function T(X). Employing this relationship, it is then proved (Map D) that 
there is a one-to-one relationship between the wave function T(X) and the corresponding 
nondegenerate ground state density p(r). Thus, knowledge of p(r) determines n(r) to within 
a constant. Since for a hxed electron number N, the kinetic T and electron-interaction 
potential W energy operators are known, so is the system Hamiltonian. Solution of the 
corresponding Schrodinger equation then leads to the wave functions T of the system. It 
is the one-to-one relationship between the external potential and the gauge invariant density 
that defines the latter as a basic variable. As the wave function T, and hence energy Efip] 
are functionals of the density p(r), the variational Euler equation for the density with fixed 
n(r) follows subject to the constraint of known electron number N (see Table 1). (The lowest 


nondegenerate 


I, llO| excited state density p®(r) of a given symmetry different from that of 
the ground state is also a basic variable.) 

In the added presence of an external magnetostatic held B(r) = V x A(r), where A(r) 
is the vector potential, the Hamiltonian when the interaction of the held is only with the 
orbital angular momentum is 


H = -Y' 

2 ^ 
k '- 


Pk + -A(rA,) 
c 


+ W + V. 


( 1 ) 


When the interaction of the magnetic held is with both the orbital and spin angular mo¬ 
mentum, the Hamiltonian is 

2 1 


^ = - V 

2 ^ 

k 


Pk + -A(rfc) 
c 


+ fT + f + -VB(rfc) 




( 2 ) 
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where s is the electron spin angular momentum vector. In deriving the Hamiltonians of Eqs. 
(1) and (2), we have hewed to the philosophy [n| that the only ‘fundamental’ interactions 
are those that can be generated by the substitution p —t p + ^A. (This then dehnes the 
physical momentum operator in the presence of a magnetic held, and thereby the physical 
current density j(r).) In non-relativistic quantum mechanics, the Hamiltonian of Eq. (2) is 
derived llj] by Schrodinger-Pauli theory for spin ^ particles via the kinetic energy operator 
icT ■ (p -(- A)cr ■ (p -f A), where cr is the Pauli matrix, and s = ^cr. The spin magnetic 
moment generated in this way has the correct gyromagnetic ratio g = 2. 

It would appear that one could prove a one-to-one relationship between the gauge invari¬ 
ant properties {p(r), j(r)} and the external potentials {n(r), A(r)} along the lines of the HK 
proof. However, no such proof is possible as the relationship between the external potentials 
^r), A(r)} and the non-degenerate ground state wave function d'(X) can be many-to-one 


12l | and even infinite-to-one |^. Hence, in these cases, there is no equivalent of Map C, 


and therefore the original HK path is not possible. The proof that {p(r), j(r)} are the basic 
variables must then differ from the original HK proof. Furthermore, the proof must account 
for the many-to-one rela 


In the literature 


donship between {f(r), A(r)} and d'(X). 

121, ll^ , the proofs of what properties constitute the basic variables are 
not rigorous in the HK sense of the one-to-one relationship between the basic variables and 
the external potentials {v. A}. Further, they do not account for the many-to-one relationship 
between {n. A} and T. Additionally, the system angular momentum is not considered. The 
choice of the basic variables is arrived at solely on the basis of a Map D-type proof between 
these assumed properties and the nondegenerate ground state T, thereby the claim that 
T is a functional of these properties. In these proofs, the existence of a bijective Map 
C is implicitly assumed, 15l. Il6| (see also last reference of 12). For example, in spin- 
DFT 2, 2 , 14| for which the Hamiltonian is that of Eq. (2) with the held component 
of the momentum absent, the basic variables are assumed to be {p(r), m(r)}, where m(r) 
is the magnetization density. In current-DFT 1^, corresponding to the Hamiltonian of 
Eq. (1), the basic variables are assumed to be p(r) and the gauge variant paramagnetic 
current density jp(r). For the Hamiltonian of Eq. (2), the basic variables are assumed 
to be {p(r),m(r),jp(r)} or {p(r), m(r), jp(r), jp,n,(r)} where jp,m(r) are the gauge variant 
paramagnetic currents of each component of the magnetization density. Subsequently, a 
Map D proof is provided. Additionally, with the basic variables now assumed known, a 
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Percus-Levy-Lieb (PLL)-type proof [^, ^ can then be formulated [l9(]. More recently, we 


gaye a derivation 


15 


20l | which purported to prove that j were the basic variables 


but the proof was in error 2l|. Subseque ntly , we proved 22| for the Hamiltonian of Eq. (1) 
that for the signihcant subset of systems jl3l. |23| for which the ground state waye function 'h 
is real, the basic variables are {p(r), j(r)}. Our proof of bijectivity between {p(r), j(r)} and 
{n(r), A(r)} explicitly accounts for the many-to-one {n(r), A(r)} to d' relationship. This 
proof then constitutes a special case of the more general proof for T complex presented in 
this work. 

Here we extend the HK theorems to systems of electrons in external electrostatic 
S{r) = — Vn(r) and magnetostatic B(r) = V x A(r) helds with known electron num¬ 
ber N and angular momentum J. The proofs are for a uniform magnetostatic held, and 
for Hamiltonians in which the interaction of the magnetic held is (i) solely with the orbital 
angular momentum (J = L), and (ii) with both the orbital and spin angular momentum 
(J = [L and S]). We prove, in the rigorous HK sense, that for fixed N and J the ba¬ 
sic variables are the gauge invariant nondegenerate ground state density p(r) and physical 
current density j(r). In other words, knowledge of {p(r),j(r)} determines the potentials 
{n(r), A(r)} to within a constant and the gradient of a scalar function, respectively. Hence, 
with the Hamiltonians known, solution of the respective Schrodinger and Schrodinger-Pauli 
equations lead to the wave functions of each system. The pro of is for (n, A)-representable 
{p(r),j(r)}. The extension to the Percus-Leyy-Lieb (PEL) [l7|, ll8| constrained-search path 
for A-representable and degenerate states readily follows. As the wave function T is a 
functional of {p(r), j(r)}, theories of electronic structure based on {p(r), j(r)} as the basic 
variables can then be formulated. 


II. PROOF OF GENERALIZED HOHENBERG-KOHN THEOREMS 

To accentuate the role of the density p(r) and physical current density j(r), we rewrite 
the Hamiltonians of Eqs. (1) and (2) in terms of operators representative of these gauge 
invariant properties. The Hamiltonians can then be written, respectively, as 

H = f + W + VA, (3) 


5 














H = f + W + VA-j m(r) ■ B(r)cir, 
where the total external potential operator Va is 

Va = '^ + ^ y j(r) ■ A(r)cir - ^ j p{r)A‘^{r)dr, 
and the corresponding energy expectations E =< \['(X)|hf|'l/(X) > as 


E — T -\- Eee + Va, 


E = T + Eee + VA- J m(r) ■ B(r)(ir, 


where the total external potential energy Va is 

Va= < d'(X)|14|d'(X) = J p{r)v{r)dr + ^ J j(r) ■ A(r)dr - J p{r)A‘^{r)dr, (8) 

and where T and E^^ are the kinetic and electron-interaction energy expectations. In the 
above equations, the physical current density j(r) is dehned in terms of the physical mo¬ 
mentum operator (p -|- ^A) as 

j(r) = N^^J vl/*(ra,X'^-i)(^p + ^A(r)^d/(ra,X^-i)dX^-\ (9) 

or equivalently as the expectation of the current density operator j(r); 


where 


j(r) =<vI/(X)|j(r)|vl/(X)> 


j(r) = jp(r) +jd(r), 


with the paramagnetic jp(r) and diamagnetic jd(r) operator components dehned, respec¬ 
tively, as 

jp(r) = -r) + 6{rk- r)pfc], (12) 


jrf(r) =p(r)A(r)/c, 


with the density operator p(r) being 
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The magnetization density m(r) is the expectation 


m(r) =< \l/(X)|rii(r)|T(X) >, 


( 15 ) 


with the local magnetization density operator m(r) dehned as 

m(r) =--^Sfc(5(rfc - r). (16) 

(The current density operator j(r) can also be dehned in terms of the Hamiltonian H as 
j(r) = cdH/dA. This conhrms that for both the Hamiltonians of Eqs. (3) and (4), the 
physical current density is the orbital current density.) 

We hrst present the proof of bijectivity between {p(r), j(r)} and {u(r), A(r)} for spinless 
electrons corresponding to the Hamiltonian of Eq. (1) or (3) for hxed electron number N 
and angular momentum L. The proof is by reductio ad absurdum. Let us consider two 
different physical systems {u, A} and {u', A'} that generate different nondegenerate ground 
state wave functions T and T'. We assume the gauges of the unprimed and primed systems 
to be the same. Let us further assume that these systems lead to the same nondegenerate 
ground state {p(r), j(r)}. We prove this cannot be the case. From the variational principle 
for the energy for a nondegenerate ground state, one obtains the inequality 


E= <T|iL|T> < <T'|.H|T'> 


(17) 


Now 


< T'|/7|T' >=< T'|T + W + V' + - j'(r) ■ A'(r)dr 


— J p(r)A'2(r)dr|T' > + < - H'|T' > 

- < T'l /[j(r) ■ A(r) - j'(r) ■ A'(r)]dr|vl/' > 


2c2 


< T'l / p(r)[A2(r) - A'\r)]dr\^' > . 


(18) 


Employing the above assumptions, and following the same steps as in 22|, one obtains the 
inequality 

E + E' < E + E’ + J [jp(r) - jp(r)] ■ [A(r) - A'(r)] dr, (19) 

where E' =< >. 
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As the majority of the experimental and consequent theoretical work is performed for 
uniform magnetic helds, our proof too is for such helds. 

Consider next the third term on the right hand side of Eq. (19). With B(r) = i?i^, 
B'(r) = B'lz-, and the symmetric gauge A(r) = iB x r, A'(r) = iB' x r, this term may be 
written as 

^=^AB-/ rx jp-jp(r) dr, (20) 

where AB = {B — B')iz. First consider the integral 


Ji = y r X jp(r)dr 


(21) 

k 

(X) [r X Vr,<5(r - r^) + <5(r - Vk)r x Vr,] ^(X). 

(22) 

Next consider the second integral of Ji of Eq. (22): 


- \) 

f dXT*(X)(^rfcXpfc)T(X) 

h 

(23) 

^0 

K 

f dXT*(X)5^LfcT(X) = ^L, 

k 

(24) 


where Lfc = x is the canonical orbital angular momentum operator, with p the canonical 
momentum operator (p = pkinetic + Pfield = mv + ^A), and L the total canonical orbital 
angular momentum dehned by Eq. (24). Note that the canonical angular momentum is 
gauge variant. 

The hrst integral of A of Eq. (22) is 

Ill = [ c?r^*(X)6«^.,^(V(r-rfc)T(X)). (25) 

k ^ v/ k'y 

On integrating the inner integral by parts and dropping the surface term, one obtains 

Ill = j dX[-eafiy j dr^^-^r;3(5(r - rfc)T(X)] (26) 

= (27) 

On integrating by parts again, one obtains 

/ o 

<'X>1-*(X) —(rt^>l.(X)) (28) 

K 

= - W f <iX'I-*(X) (r, X V,J'I-(X) = 1l 
k d 
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(29) 



Hence, the integral I of Eq. (20) is 


2 ^ ’ 


(30) 


If one imposes the condition that the total canonical orbital angular momentum is fixed so 
that L = L', then the integral I vanishes so that the third term on the right hand side of 
Eq. (19) vanishes. 

For states with fixed orbital angular momentum L, Eq. (19) then reduces to the contra¬ 
diction 

E + E' < E + E'. (31) 

What this means is that the original assumption that d' and \h' differ is erroneous, and that 
there can exist a {v. A} and a {n'. A'} with the same nondegenerate ground state wave 
function. The fact that T = T' means that p(r)|,i, = p'(r)|,i,/. However, the corresponding 
physical current densities are not the same: j(r)|,i, ^ j^r)!^', because jd(r)|^ ^ if one 
hews with the original assumption that A(r) is different from A'(r). This proves that the 
assumption that there exists a different {n'. A'} (with the same N and L) that leads to the 
same {p,j} as that due to {v. A} is incorrect. This step takes into account the fact that 
there could exist many {n. A} that lead to the same nondegenerate ground state T. Hence, 
there exists only one {v. A} for hxed N and L that leads to a nondegenerate ground state 
{p, j}. The one-to-one relationship between {p, j} and {n. A} is therefore proved for the case 
when the interaction of the magnetic held is solely with the orbital angular momentum. 

With {p(r),j(r)} as the basic variables, the wave function T is a functional of these 
jroperties. By a density and physical current density preserving unitary transformation 
^,ll5|,l2^ it can be shown that the wave function must also be a functional of a gauge function 
a (R). This ensures that the wave function when written as a functional: T = \h[p,j,Q;] is 
gauge variant. However, as the physical system remains unchanged for different gauge 
functions, the choice of vanishing gauge function is valid. 

As the ground state energy is a functional of the basic variables: E = E.„^a[p, j], a varia¬ 
tional principle for E,^^x[pi'f\ exists for arbitrary variations of (v, A)-representable densities 
{p(r),j(r)}. The corresponding Euler equations for p(r) and j(r) follow, and these must be 
solved self-consistently with the constraints / p(r)(ir = N, J r x (j(r) - ip(r)A(r)dr = L 
and V ■ j(r) = 0. Implicit in this variational principle, as in all such energy variational 
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principles, is that the external potentials remain fixed throughout the variation. (See Table 

I.) 

We next consider electrons with spin corresponding to the Hamiltonian of Eq. (2) or (4). 
In this case, with the same assumptions made regarding the two different physical systems 
{n. A; -0} and {f'. A'; leading to the same {p(r), j(r) as before, the inequality of Eq. (19) 
is replaced by 

E + E'< E + E' + j [j;(r) - jp(r)] • [A(r) - A'(r)]o!r 

— I [m'(r) — m(r)] • [B(r) — B'(r)] dr. (32) 


The third term on the right hand side vanishes if one imposes the constraint that the orbital 
angular momentum L of the unprimed and primed systems are the same. Hence, next 
consider the last term of Eq. (32). For a uniform magnetic held with B(r) = Bi^ and 
B'(r) = Biz, we have 


where 



J m(r) ■ B(r)dr = B 


mz{r)dr, 


mz{r) = -^[pa{r) - P0{r)], 

with pa{r), pgi^r) being the spin-up and spin-down spin densities, 
inequality is then 


(33) 

(34) 

The last term of the 


j [m'(r) - m(r)] ■ AB(r)dr = j [{^^(r) - ^(^(r)} - {p«(r) - p^(r)}] dr, (35) 

with AB = B — B'. If the ^-component of the total spin angular momentum Sz for the 
unprimed and primed systems are the same, the corresponding spin densities are the same. 
The last term of Eq. (35) thus vanishes leading once again to the contradiction E + E' < 
E + E'. More generally, the magnetization densities m(r) and m'(r) are the same if the 
total spin angular momentum S are the same. Hence, once again, the bijective relationship 
between the nondegenerate ground state densities {p(r), j(r)} and the potentials {t'(r), A(r)} 
is proved provided one imposes the constraint that the total orbital L and spin S angular 
momentum are hxed. 

This may be seen in a different manner by accentuating the role of the spin angular 
momentum. With the 2 ;-component of the total spin S being Sz = Yhk the density 
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m^(r) may be written as 


(36) 


mAr) = - - ^ S'^ 7 (rcr, ra), 


cN 


with 7 (xx') = N J 'h*(rcr, ^)'h(rV', ^)(iX^ the density matrix. Since in the 

primed system, the spin vectors are different, i.e. some s(,, we have 

J [m'(r) — m(r)] ■ AB(r)dr = AB J — mz{r)~\dr (37) 

= [>5'7'(ra,rcT) - ^;, 7 (ra,ra)]dr. (38) 

<7 

Employing the original assumption that the diagonal matrix elements 7 (r(T, ra) of the density 
matrix 7 (xx') are the same for the unprimed and primed systems we have the right hand 
side of Eq. (38) to be 

^ y Sz]-r{ra, ra) = 0 (39) 

(7 

provided S'' = S'^. 

In the above proofs for the Hamiltonians of Eqs (3) and (4), the definition of the current 
density j(r) employed is that of Eq. (10). However, for finite systems, the Hamiltonian of 
Eq. (4) can also be written as [25 1 


H = f + W + V + ^ y jp(r) ■ A(r)dr + ^ /p(r)H^(r)dr 


1 


jm(r) ■ A(r)cir, 


(40) 


where the magnetization current density operator jm(r) is defined as 


jm(r) = -cV X m(r) 


Hence the physical current density j(r) may also be defined as 25 1 


(41) 


j(r) = = jp(r)+jd(r)+j™(r), (42) 

the sum of the paramagnetic, diamagnetic, and magnetization current densities. Even for 
this definition of the physical current density j(r), the proof of bijectivity between {p, j} and 
{n. A} is valid provided the angular momentum L and S are hxed. For spin-compensated 
systems, the magnetization current density jm(r) vanishes. 
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III. CONCLUDING REMARKS 


In conclusion, we have generalized the HK theorems to the added presence of a uniform 
magnetic held. We have considered the cases of the interaction of the magnetic held with 
the orbital angular momentum as well as when the interaction is with both the orbital and 
spin angular momentum. In this work we have proved a one-to-one relationship between the 
external potentials {n(r), A(r)} and the nondegenerate ground state densities {p(r),j(r)}. 
The proof dihers from that of the original HK theorem, and explicitly accounts for the many- 
to-one relationship between the potentials {n(r), A(r)} and the nondegenerate ground state 
wave function T. To account for the presence of the magnetic held, which constitutes an 
added degree of freedom, one must then impose a further constraint beyond that of hxed 
electron number N as in the original HK theorems. For the Hamiltonian corresponding to 
spinless electrons, the added constraint is that of hxed canonical orbital angular momentum 
L. For that corresponding to electrons with spin, the constraints imposed are those of hxed 
canonical orbital L and spin S angular momentum. (The gauge employed for the canonical 
angular momentum L can be chosen to be the same as that employed for the Hamiltonian.) 
It is the further constraint on the angular momentum that makes a rigorous HK-type proof 
of bijectivity between the gauge invariant basic variables and the external scalar and vector 
potentials possible. Additionally, the HK-type proofs are possible because the Hamiltonians 
considered are rigorously derived from the tenets of nonrelativistic quantum mechanics. 

With the knowledge that the basic variables are {p(r),j(r)}, a variational principle 
for the energy functional for arbitrary variations of (n, A)-representable densi¬ 

ties {p(r), j(r)} is then developed for each Hamiltonian considered. The constraints on the 
corresponding Euler equations are those of hxed electron number and angular momentum, 
and the satisfaction of the equation of continuity. 

Again, knowing what the basic variables are, it is possible to map the interacting system 
dehned by the Hamiltonians of Eqs. (1) and (2) to one of noninteracting fermions with the 
same p(r), j(r), and J. Such a mapping has been derived within QDFT 2^. The theory 
has been applied to map an interacting system of two electrons in a magnetic held 
and a harmonic trap n(r) = for which the ground state wave function is \['(ri,r 2 ) = 

C{1 -|- ri2)e“5hi’'2)^ where ri 2 = |ri — r 2 | and C'^ = l/7r^(3 -|- \/^), to one of noninteracting 
fermions with the same {p(r),j(r)}. This example corresponds to the special case of zero 
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angular momentum. However, the QD 


forward. For other recent work see 


27 


T mapping for hnite angular momentum is straight 


28| . The conclusions in the latter are based on the 


assumption of existence of a HK theorem but one without the requirement of the constraint 
on the angular momentum. 
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Theory 

Hohenberg-Kohn DFT 

Generalized HK DFT 

Parameters characterizing 

ground state 

Electron Number N 

Electron Number N 

Angular momentum L 

Relationship between 

potentials and wave function 

One-to-one between 

n(r) and T 

Many-to-one between 

{n(r), A(r)} and 'k 

Properties characterizing 

ground state 

Electron density p{r) 

Electron density p(r) 

Physical current density j(r) 

Bijectivity theorem 

For fixed N 

p{r) -H- n(r) 

For fixed N and L 

{P(r),j(r)} O {n(r), A(r)} 

Wave function 

and Energy functionals 

T = ^[p,a] 

For fixed v : E = Ey[p\ 

T = ^[p,i,Oi] 

For fixed {n, A} : E = 

Euler equations 

and constraints 

Variational principle for 

fixed V and known N: 

&Ev\p] _ n 

5p - U 
f p(r)dr = N 

Variational principle for 

fixed {n, A} and known A, L: 

A [Pij] (-1 c\ 

6p - ^ 6} - ^ 

J P 

f p(r)dr = N 

frx (j(r) - ip(r)A(r))dr = L 

V • j(r) = 0 


TABLE I: Comparison of Hohenberg-Kohn and Generalized Hohenberg-Kohn theories. 
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